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1. TNTRODUCTION 


A. -DESCRIPTION TOPS THE SEROR EECN 

It is very common for users to access the same channel 
using the Time Division Multiple Access (TDMA) or Frequency 
Division Multiple Access (EDMA) techniques. These two 
methods provide multiple accessing, however there are disad- 
vantages. For example reassigning of channel is difficult 
in FDMA. In TDMA guard times and headers reduce throughput. 
TDMA also requires central timing and synchronization 
control. In addition, these two methods are not suitable for 
handling interactive traffic such as that carried by 
computer networks. For this type of (traffic pr l ep nt o K D 
peak load to average load is rather high. Users generating 
such traffic are often called bursty users. 

An approach suitable for burstv users is to access the 
channel randomlv. This is known RANDOM MULTIPLE ACCESSING 
(RMA or RA). In a random multiple access svstem each user 
has equal opportunitv to access the channel when he has a 
message to send. Naturallv, there is no coordination between 
users and collision is possible. Two or more messages trving 
to access the channel at the same time will lead to overlap- 
ping of messages which is called a collision. Overlapping 
of messages results in partial or complete destruction of 
each message and is not desirable. Collided messages need 
to be retransmitted until thev are successfullv received by 
their receivers. In the next section, we shall provide a 


brief historical review of random access techniques. 


Deeds CAL TOR OF THE RANDOM ACCESS TECHNIQUES 
l. Pure ALOHA 

In September 1968, the University of Hawaii 
[Ref. 1], began a research program to investigate alterna- 
tives to the use of conventional wire communications for 
computer-computer and console-computer links. The work was 
conducted by Norman Abramson and resulted in a system which 
is now known as the ALOHA system. 

ALOHA was a system for remote stations on various 
islands of Hawaii to contact a central computer via a common 
radio broadcast channel. Stations use a very simple protocol 
for accessing the channel. This protocol is known as pure 
ALOHA. 

The fundamental idea of pure ALOHA is very simple. 
Whenever a station or a user has a message it is transmitted 
immediately. No coordination is required for the user. 
Messages are first sent to the central station and then 
repeated back to all users in broadcast mode. In other 
words, the central station just acts as a repeater. Two 
possible outcomes may result from the transmission, success 
or collision. A success means the message does not overlap 
with other messages in time and will be successfully 
received by the receiver. Collision leads to destruction of 
all the messages involved in the collision. Therefore, after 
transmission of the messages the sender should monitor the 
feedback channel. If collision is detected, which can be 
done through parity checking, retransmission of the collided 
messages is required. In order to prevent unending colli- 
sion, retransmission should be done in the following manner. 
Once a retransmission is necessary, a randomized delay is 
usually introduced before the retransmission actually takes 
place. Figure 1.1 shows of the execution of this protocol. 
Mom packets C and D collide during the fisst 
attempt. 





Time 


Fiġure mm Pure Aloha 


ALOHA protocol is its simplicitv. A disadvantage is that it 
can only provide very low stable throughput. In this case, 
throughput is defined to be the average number of packets 
that can be successfully transmitted during one packet 
transmission time. A notation that has been generally 
accepted kar Ehroughpae in. Clearly O < = >) lira 
evaluating S, we usuallv assume an infinite number of users 
so that the channel input can be modeled as a Poisson 
process with mean arrival rate G packets/ 7 sec. Here 7 is 
the packet transmission time in seconds. Channel input 
includes newly generated and retransmitted packets. It can 


be proved that, 
> = G exp(- 25) 


The maximum of S occurs at “G=G.5 Which se me 5 Mika 2 
throughput 0.18. A system which can only achieve a maximum 
throughput if 0.18 is not efficient At Alt one rèéasoòat gon 
this low efficiency is due to a 27 second vulnerabmilev 
period of a transmitted packet. This situation is depicted 


in Figure 1.27 


10 


Vulnerable period of packet * 





Araure 1.2 Vulnerability Period of the Shaded Packet 


2. Slotted ALOHA 
One way to improve the performance of the pure ALOHA 


system is to reduce the vulnerability period of the packet. 


Ime 1972; Roberts (Ref. 2], introduced the slotted ALOHA 
protocol for this purpose. The idea is demonstrated in 
BRE 6ê 6 2 Im Which the channel time axis is sliced into 
T-second intervals called slots. Retransmission of data 


packets can take place only at the beginning of time slots 
as shown in Figure 1.3. Clearly, the vulnerability period 
BBS Packet 15 now T Seconds, which is one half of that of 
the pure ALOHA protocol. However the price we pay here is 
the requirement for strict synchronization among all the 
participating users. Once again, assuming an infinite number 


of users it can be shown that, 
S = G » expl(-G) 


The maximum throughput becomes 0.36 which is twice 
Mao ene pure ALOHA protocol. Figure 1.4 shows the 
(AMOU HDUESEnaracteristics of pure and slotted ALOSA, 

Although the system capacity, i.e., the maximum 


achieveable throughput is doubled by slotted ALOHA, th 


ປ 





ະປະສ S 10 io 20 


Eigumeė Lo Slotted Aloha Channel 


system is still not stable. If the mean arrival at G is 
pushed over 1.0 then there will be a drastic increase in the 
coll usina tes This will eventually lead to a deadlock 


situation when the system throughput is reduced to zero. 


0.4 
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Ala AC củ 


6.0 8.0 10.0 


Figure 1.4 Aloha Svstems 
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3. Collision Resolution Algorithms 
Since the introduction of the ALOHA protocol by 


Abramson, tremendous efforts have been invested DY 
researchers trying to enchance the performance of a random 
access system. Numerous papers have appeared as a result of 
these efforts. These recently have been centered around the 
collision resolution algorithm (CRA). A collision resolu- 
tion algorithm is a protocol which is designed to reschedule 
the retransmission of those users involved in a collision so 
that each user eventually gets through. Usually before a 
collision is completely resolved transmission of new packet 
should not be attempted. The efficiency of a CRA is usually 
measured by the average number of time slots required to 
resolve a collision. This average number of time slots is 
usually called a collision resolution interval (CRI). Based 
on measurement one can usually derive the maximum achieve- 
able throughput provided by the protocol. Almost all CRAs 
use a slotted system, 1.e., the channel time is slotted. 
Historically, the first CRA was proposed by 
Capetanakis [Ref. 3]. A similar approach was proposed by 
Hayes [Ref. 2]. Capetanakis' collision resolution algorithm 


(CCRA) was later improved by Massey [Ref. 4], by Gallager 


[Ref. 2] and by Humblet [Ref. 2]. An energy detector was 
also employed CO improve the performance of the 
Capetanakis's protocol in IRef. 3). A thorough review of 


CRAs is provided in reference 2. 


We use the example sketched in Figure 1.5 to illus- 
trate Capetanakis's idea. Here the time axis is organized 
into 2 r7-slot frames. At the termination of the previous 
collision resolution interval all of the users are waiting 
for transmission of their packets. Each one randomly picks a 
slot in the next frame and transmits his packet. In Figure 
1.5 we have a total of five users trying to access the 


channel. They are represented by A, B, C, G and H. Users A, 


15 
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Figure 125 Binary Tree Algorithm 


B and C select the left slot, while G and H use the right 
slot. Of course collisiom 'occurs' 1: 125 E lo Esa Slot status 


information available to users is ternary, i.e.,empty, 
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EM ele MOD collision. Since a collision is detected in slot 
1, we then proceed to resolve the collision. We let users 
involved in the collision of slot 1, independently and 
randomly select slot 3 or 4 and retransmit. By this we 
F i onu ridae A, B and C into two groups. On our example 
ce select” slot 3, and € uses slot 4. Since C is now a 
success, i.e., then proceed to resolve the collision between 
A and B. Collision between G and H can be resolved only when 
A, B, and C have been successfully transmitted. A tree 
representation of Capetanakis protocol is given in Figure 
mea). 


Tree Algorithm : Each leaf on the tree represents a source 
and each source has three bits binary address in this 
example. It might be 8 sources in the tree. There are five 
active sources in the Figure 5(a). The root of the tree 
defines all collision resolution interval, First upper 
branch leaves go to the left slot of the first group (A, B, 
C in the slot 1). Lower branch leaves go to the right slot 
of the first pair (G, H in the slot 2).The system shows the 
retransmission of the A, B, C in the second slot pair, C is 
alone and means success. After retransmission packet A and 
ວໃ slots 5 and 6. Up here CCRA finished the upper 
branch collisions. Next step is to retransmit G and H by 
using the same way. 

Capetanakis provides two versions of his protocol, 
Static and Dynamic. The maximum achieveable throughputs are 
0.346 and 0.430, for static and dynamic CCRA respectively. 

Later, Massey improved Capetanakis's algorithm by 
observing that if a collision is followed by an empty slot, 
one slot can be saved by repeating the random retransmis- 
sion. This modification allowed a throughput of 0.46. 

Gallager introduced a different conflict resolution 
algorithm with guaranteed stability and had first come first 


served characteristics. He selects an initial arrival 


Hi 


imterval. ` P. IL performs subsequent subdivisions of 
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whenever collision occurs, and it reinitializes whenever a 
collision occurs Within the first ENBAIVISTONS e m mecen 
obtained a maximum throughput of 0.4872. Humblet improved 
Gallager's. algorithm “and increased -nc m: ou pu Mo 
0. 48/75: 

In 1982, Georgiadis and Papantoni-Kazakos [Ref. 5] 
proposed a.collision, resolution" aličjo muhi ca) r 
CRAI employs a device called an energy detector to tell the 
number of the users involved in a collision. Based on this 
information, they suggest an optimal way of dividing users 
into two groups once a collision involving these users is 
detected. Their protocol offers a maximum channel throughput 
of 0.53. The difficulty with this protoco Mi - u E 


assumes the availability of an infinite energy detector. 


C. TRE*SCOPE (OM slat tomo o 

In this thesis we examine two new protocols which are 
not only easy to implement but also offer satisfactory 
performance. In Chapter II, we study the MSCRA protocol 
where each, time we open a number of time slots and let 
users involved in the collision randomly choose one slot and 
retransmit. The number of slots opened equals the number of 
users involved in the collision. Therefore, we have to 
assume the availability of an infinite energy detector. 
Chapter III gives an adaptive version of this protocol. 


Finally, the conclusion is presented in Chapter IV. 
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IT. MULTISLOT COLLISTON RESOLUTION ALGORITEM 


A. INTRODUCTION 

Suppose we have an energv detector of infinite capacitv 
so that when a collision occurs, we know exactly how many 
users are involved in the collision. After the detection of 
a collision involving n users, we immediately open the next 
ອ ots. for collision resolution. We let each of the 
collided users pick one of these m slots randomly for 
Mere nemi sskon, and ' m " should be selected to increase the 
average number of successful transmissions in the selected m 
slots and to maximize throughput. Towards this goal, we 
shall now prove that the optimum m happens to be n (n is the 
number of the collided users). 


Suppose we wish to distribute n packets into m slots. 


X = The random variable representing the number of the 
slots which contain exactly one packet when 
dlstributing'n packets imto m slots. Physically X 


represents the number of successes. 


Next let E[X] denote the expectation of X. Then we have the 


following special cases. 


Case l-: me il mel ; ELX] = 1 
Case 2 : n=2 m=2 ; It is impossible to have only one 
success in this case. Distributing two packets into two 


slots, which is equivalent to throwing two balls randomly 
into two boxes, always results in either two successes or no 


success at all. In other words, 


PE — 240173) + 0*(1/2) = 1 
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Case 3 : n=2 m=3 


m m-l 3 2 
E[x] = 2 / mé =2 ກ GO 
າ 1 1 1 
Case 4 : n=2 IA 2 
m m- 1 
= A = ši 2 
ELE) = 2 mi = 2m(m-1)/m 

d 1 


The denominator mé represents the number of ways to throw 2 
balls into m boxes. The numerator tells the number of ways 
of select two distinct boxes from m boxes so that each box 
containing exactly one ball. RE 


of m. To find the m which maximizes E[X] we perform 





d E d 2(m-1) 2m-2m(2(m-1)) 
— < — = = 
dm m dm m? m 


which leads to m=2 
Case 5 


n=3 m=3 In this case we never can have exactly two 
successes when distributing 3 packets into 3 slots. If the 
number of successes is two, the third packet must also be a 
success. Therefore we can have either one or three 


successes. Thus, 


E[X] = 1 ເຈ 27: DASS 


iko 


Case 6 


n=3 m>3 


m /A fm-1\ fn-1 m /m-N /m-2 3m(m-1)l 
E[X] =|1 +3 / mě? < = — 
1/\1/\ 1 2 IA 1 1 m? 


In erdèr to find optimal m, again we carry out 





IMMA EX] a [(m-1)2] m°2(m-1)-(m-1)?m* (m-1)(3-m) 
3 dm m dm m? m? mé 
(2783 


which vields m=3 


In general, it can be easilv seen that 


nm(m-1)n 1 n(m-1)»71 


FU ແນ 


If we let m approach infinitv, the l'Hopital's Rule shows 


n(m-1)271 n(n-1)(m-1)"~4 
EE ເຢ = Lim = lim 
` 2 „nl ka“ (n-1)m2~2 
n(m-1yR~2 
DE L = n 
KV nh = e 


which is intuitively reasonable, for the optimal selection 


of n we have. 
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1 £ 3 າ (n-1)(m-1)?~2m1-nmR" ?(m-1) 8" 1 


n dm m mén 


ml (m-1)""4(-m+n) 





2n 


Which results in m=n. 
B. DESCRIPTION OE THE MSCRA 


First, we outline the fundamental assumptions to be used 


im thisjehap sen? 


de Traffic generated by network users collectively 


form a Poisson process. 
25 Channel time is slotted. 


Še An energy detector of infinite capacity assumed in 


reference 3 is available. 


4. There are no new packet transmissions until the 


collision has been resolved. 
> There is mo propagation delay, 


6. There is no channel errors. 


Based on our discussion in Section II.A the MSCRA can De 


stated as follows. After a collision involving n users the 
next n slots are opened for retransmission. Each user 
randomly selects one of the n slots and retransmits. At the 


end of the n slots if i of the users have successfully 
transmitted their packets, i.e., n-i of them are still 


involved in collisions, the process is repeated atenn 
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EmWčure 2.1 multislet'Collision Resolution Algoriehm 


In this example, the energv detector shows 5 collided 
packets, at the verv beginning and the svstem opens the next 
Emvemslots for collision resolution. During this five-slot 
period there is only one success and four packets still 
collided. The system then opens four more slots. Now A and 
B are in collision while D and E succeed. It is thus neces- 
wero open another two slots (slot ll and 12). In our 
example A and B are now successes. 

It can be seen from this example that after each retran- 
smission, the number of successes and still colliding 
packets can be assertained. The number of yet to be opened 
slots is always equal to the number of the packets remaining 
to be resolved. In this example, the length of the CRI 


(Collision Resolution Interval) happens to be 12 slots. 


Ce eoNALYSis OF THE MSCRA 

In this section we begin by finding the probability of 
weona Successes in distributing n balls into m boxes. 
ວ... ເ] ) . 


SR 


il 
jou 


if n=l and m-ce OoN ncl P(X-1) 


no oe 2 
if n=2 and m=2 then P (X=0) = ແ (2) <... > 
PIYE) O 
2 SÀ Ni l 7 =0. 
POSI = \( ( 2 
Ji il l iL 


if n=3 and m=3 then 


3 8 
1 
Peo) -( )( ) ” ER 
SE 3 
3 a 2 2 
3 2 
Pen = (JÍ Il )( ), 3 = 
1 jů l2 
ບ ບ > 0 
2 
P(x=3) =— 779" 


In resolving a collision involving n users, we define 
the system state to be the number of collided packets 
remaining to be resolved. Figure 2.2 shows the state tran- 
sition diagram in resolving n collided packets. The tran- 
sition from state n to state i is labeled by a transition 
n,n-i' THUS Pc n-i 


of having n-i success in an n slot period. 


probali a represents the probability 
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ຈຄາ. 2 . 2 State Transition Diagram of An n slot Period 


An important performance measure of the protocol is L, 
the expected number of slots required to resolve a collision 


involving n users. From Figure 2.2. 


Ly = Pn, on ħi 2n,12n-1 +...+ P 


EN d Ph 


mn „no (2.2) 


Where Lo = Li z 0 


Since it is not possible to have exactly n-l successes, we 


have 
Ph,n-l = 0 
F l Tg Eqn (2.2) 
n n 
Ly = A > | Pa, ibn-i ມ, Pa obn +) k n-i 
2 — 0 1 = 1 


We obtain, 
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l 
ji = — ງ r. = te) 
n 1l. m, 1 Ta h 


Finally, the expected number o£ slots required in transmit- 


ting a total of n packets is, 
A a 
Where the vis due to the slot of the inita e e en. 


P can be obtained recursivelv as follows, 


AR E 
Pai = À ວຫາ / n 
Where 
An + The number of possible arrangements in distrib- 
4 
Uting n packets into NM g CS, which sesult in 
successes. 
n SA j l l E 
+: : ສ. n= PSY = 
a ` l 1. (n-i) - A i -A 2. 4G . f 
id A Dei a 55" “2 ດ 24 


Eqn (2.4) is valid for © < 7 =: “+ 


n n 
n = ' = ' 
Ann mi n! 
m I) 


D. NUMERICAL CALCULATIONS AND DISCUSSIONS 


For the MSCRA protocol two major computed programs were 
written. The first one provides numerical values for E of 
Eqn (2.3). The other program givesWšimulai lon results ae. 


verify me 
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ie: Numeracal Calculations of i 


ມນສ > program, wasygysittengin WATEIV using doubl 


iD 


precision and real value variable declarations. This program 


Pecursavely generates the ANS, ; of Eqn (2.4) in he 
4 


following order; 


AQ,o 
À1 1 A; 0 
À2 2 À2 1 À2 O 
À3 3 2302 A3 1 A3 0 
Ti 
4 E 
e š | | 
= | 
| > 
1 - 
1 
O l 2 3 4 | 
Figure 2.3 The Relationship Between n and i 


As demonstrated in Figure 2.3, i cannot be greater than n 
(the number of successes cannot be greater than the number 
ວ >. n After calculation of the ເ we can obtain L.. 
easily. A function subprogram was written to calculate the 


FACTORIAL. The Watfiv program is given in Appendix A. 
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2. demmu ila tion 
The second program gives simulation results bv using 
GGUD (a random number generator) from the IMSL routine 
librarv. GGUD is used to control the random retransmission 
of the packets. The flow diagram of the simulation is 
illustrated in Figure 2.4 . A programWvlisti nemi - Ms rno acded 
in Appendix B. 
3. Examples and Discussions 
Table I gives a comparison between the ha obtained 
by analytical formula Eqn (2.3) and that obtain by sìn 
tion, We observe that the agreement between analysis and 
simulation is extremely good. Since the error is always less 
than 1%, thus verifies the correctness of the analysis done 
for the MSCRA. This comparison is also represented graphi - 


cally in Figure 275S 
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Figure 2.4 Simulation Path 
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ANALYSIS 
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III. ADAPTIVE MULTISLOT CRA 


A. INTRODUCTION 

The efficiencv of the MSCRA can be improved bv modifving 
the algorithm such that counting for the number of 
successful transmissions can be made in an adaptive manner. 
As described in the previous chapter, the MSCRA opens n 
slots when a collision involving n users is detected. It 
then waits until the end of the n slot period, counts the 
number of successful transmissions and restarts with the 
remaining packets yet to be resolved. 

In our new algorithm, upon the detection of a collision 
involving n users we still open n slots but make observation 
on only the the first b ( b< n ) slots. At the endito - 
bth slot, we count the number of successes which have 
occured in the first b slots and repeat the process. This 
new algorithm will be called Adaptive Multislot Collision 
Resolution Algorithm (AMSCRA). 

In this chapter, we shall consider the following two 
types of AMSCRA: 

(1) If i successes have occured in the first b slots, 
repeat the collision resolution procedure with m-i Users Tin 
other words, without waiting until the end of the nth slot, 
we immediately open n-i slots and proceed similarly. We 
Shall call this policy AMSCRA with collective resolution. 

(2) If j users have tried their retransmissions in the 
first b slots and among them i have succeeded, then repeat 
the process immediately with j-i users. After each of these 
j-i users have successfully transmitted, restart the process 
again with n-j users. This policy will be called AMSCRA 


with separate resolution. 
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B.  AMSCRA WITH COLLECTIVE RESOLUTION 
ນທະບ < ເເກ 1... of the Protocol 

After a collision involving n users is detected, the 
System Opens n slots. Each of colliding users randomly picke 
ແນວ and retransmits. Users participating in the collision 
resolution observe the status of the channel for the first b 
slots. Let i represent the number of successful packets in 
the first b slots. This process is repeated at the end of 
the bth slot with n-i users. If the number of the packets 
to be resolved n, becomes smaller than b, then open n slots 
and make observations over this n-slot period. This case 
reduces to the MSCRA of Chapter II. 

In the example of figure 3.1 the energy detector 
indicates that the number of packets in the first slot is 5. 
The resolution algorithm opens five slots and observes the 
first two slots (slots 2 and 3) because b=2. Unfortunately, 
no successes are observed in slots 2 and 3 so system opens 
another five slots following slot 3 and checks 2 slots again 
(slots 4 and 5). Now we see a success for message A in slot 
4. However there are still four collided packets which need 
to be resolved. The systems opens four slots, checks slots 
6 and 7 and finds two more successes for messages B and D in 
this third interval. The total number of successful packets 
is now 3. The fourth interval is now opened and messages C 
and E are also observed to be successful. Thus the entire 
collision resolution interval(CRI) is completed. ih tnis 


particular example, the CRI requires a total of 9 slots. 


2. Analysis 
Before proceeding to derive a a few definitions are 


in order. 


gheprobability that out of n users- j of chen 


retransmit in the b slots and result in i succees. 


3 





ໄປ 


E 

B D E 
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Figure 31 An Example for AMSCRA with 


Collective Resolution b=2 


yP, > Prob[i successes/j users retransmit in the lst. b 
slots). Prob[Among n users j of them retransmit 


in thesist.b slots] 


AP. . = The mumber of the axrrangements in distributinç j 


packets into b slots, which result in i successes. 


In general we can write, 
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BR Qu [eo taba, ເ... 00 
Ban (3.2) is expressed in recursive form and is 
F iO tor OBS 1 < min(b,j)« The factor kal represents tne 
number of arrangments in selecting 1 packets from b slots 
and () repyesents the number of “mrangments possible in 
`... i packets out of j packets. The factor i! repre- 
sents the number of arrangments possible in distributing i 
meckets into i slots such that each slot is a success. In 
Ea (3.2), (b-i) represents the number of arrangments in 
distributing the remaining j-i packets into b-i slots which 


result in no successes. If i-min(b,j), then, 


k=min(b,j) 
D j 
AI = k!-U(b-j+1) ແ. 
k k 
di 20 
where U(x) = 
O otherwise 


Moe). in Egn (3.2) is used to prevent b 2 1, i=min(b,j) 


meet Occurring. Based on the definitions of vo, a and 
b 
A is we have 
A,?; mi bo j b n- j 
be | yonn" — LB off, 
Và i T າ ( 4) 
b3 3 n n 


ae mai 50 probability the system state proceeds from n to 


pee the Situation illustrated ın Figure 3.2. 
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EAI GUM 2 Transition From State neve iii 


TE ma B 
Tnus 1ni F3 
n ເວມ pa 
= za ZA 
LM = yo vn Dr STE : 
BŽ D Ya Oja tm) E 
J=0 150 
| 1 nel 
where L = © 5 
| E otherwise 


- A 
The relationship between L and WER S necesa! 


because if n-i = 1, then we need one more slot to complete 
the CRI. However L; has been EECHER 
Pana lly, the actual number of time Slots needed in 


successful transmitting nm packers ems 
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LE TW 09 as we polinuagdsout earlier im thas chapter, 


ie esult reduces to that in Chapter] lpethus 


1 n 
AN 
n + P L neb 
1-p A | n,i n-i 

P isl 

= (3:69 
n 
manb) 
l YAN 
b + AAA L 
l- V i | ດ: MA a 
3,0 1= 

150 1=O 

Where l n=] 
A | 
lè, = 
ວ 
| L. otherwise 
on Numerical Calculations and Discussions 
a. Numerical Calculation of Ly 


A COmpUutertproctammtito calculate EA was written 
in WATFIV. Double precision and real value variable declara- 
tions were used. This program creates AP. Of Son (827) ana 
Em (223) recursively. This program differs from the 
program of Chapter 2 in that we must ensure that i is never 
greater than b or j (the number of successes is limited bv 
the number of slots and the number of users). The Watfiv 
program of AMSCRA with collective resolution is given in 


Appendix A. 


YO 


b. «Simulation 
A simulation program was written in FORTRAN. 
GGUD was used as a random number generator. The Fortran 
program is provided in Appendix B. A flow diagram of the 


simulation is illustrated inc 


Read n, b 


Open n slots, call 


random number generator 


n times 





Find the number of packets appearing 


in each slot, observe the first b slots, 


count the number of successes i, 
L =L +b 
n n 
n = n iu 
no 
| stop 
| 
yes 
b < n 
no 


Figurer m: Simulation Flow Diagram for AMSCRA 


with Collective Resolution 


So 


C. AMO O S and Discussion 

Table II on p. 47 gives a comparison between the 
La obtained by analvtical formula and those obtained by 
Simulation. We observe that the agreement between analysis 
and simulation is extremely close. This verifies the 
BOMECEENEeSss of analysis done fomathis protocol. Gompazsison 
SEL Se done graphically in Figures 3.4, 3.5, 3.6 an@ 3.7. 
Moeemesrigure 3.4 we observe that b = 1 offers the Best 
performance. Comparison between MSCRA and AMSCRA will be 


made at the end of this chapter. 
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Figure 3.4 AMSCRA with Collective Resolution b=1 
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Figure 3.5 AMSCRA with Collective Resolution b=10 
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Figura 380 AMSCRA with Collective Resolution b=20 
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Figure 3.7 AMSCRA with Collective Resolution b=30 


B'SMWANSCRA WITH SEPARATE RESOLUTION 
IN DEScriprion 

After a collision involving n users is detected the 
systems opens n slots. Again each user randomly picks one 
slot and retransmits. All the participating users observe 
ste channel until the end of the first B slots j out of n 
W ess “have tried in the first b slots resulting in i 
successes. Then we repeat the process twice with j-i users 
and the other with n-j users. In other words, at the end of 
the bth slot we immediately open j-i slots to resolve the 
collision among the remaining j-i users. After each of these 
j-i users has successfullv transmitted. We open n-j slots to 
resolve the collision among the remaining n-j users. In any 
case, if n « b then we only open n slots and make observa- 
Bon for these n slots. 

curse. 3.9, the collision in slot 1 involves" pise 
users. The system therefore opens five slots for collisior 


resolution. Each of these five collided users selects one 
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Figure 3.8 AMSCRA with Collective Resolution 
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slot for retransmission randomlv. In this example A and 3 
select slot 2, D and E select slot 6 and user C selects slot 
me divide the collision reselution into two parts. the 
mẽ... part covers the first b slots (slot 2 and 3), ae 
second part covers the remaining slots. The system stores 
the users which appear in the second part (user C in slot 
4). In the first b slots we observe no successes and four 
collided packets (j=4); therefore the system opens é more 
slots, (4, 5, 6 and 7). There are now two users, A and 8, 
involved in a collision the first two slots of the 4 newly 
opened slots. The algorithm then opens two more slots (slots 
6 and 7) to resolve the collision between A and B. Users A 
and B now successfully transmit their packets in slots 6 and 
7. After this step, the system opens two slots (8 and 9) for 


the previous stored packets D and E. Both D and E now 


succeed. Finally the system opens one more slot for packet 
om 
2. Analysis 


After the observation of the first b slots, the 
Groupes of n collided users will De subdivided into two sepa- 
ae Groups withssiżes j-i and n-j, j represents the number 
Sa Users Which retransmited in the first b slots. This situ- 
on is now demonstrated in Figure 3.10 . ` ¡ has been 
mac co in Section III.B.2. 


E figure 2.10, if D 2 n we have, 


= Sp b A ZN 
ທ > ວ ເເ s (3.7) 


j=0 1=0 
l ir “2-1 
AN 
where Los 
E otherwise 
j ແ EH Sec 


In general 
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AMSCRA with Separated Resolution 


Figure 3.9 





Bigure 3.10 State hrainsitiom Diacggam for 


AMSCRA with Separated Resolution 


n 
1 2 
n + P "` « < 
1=P ) | Ay. HA n < zh 
n, 0 i=0 
L = 
n 
(35) 
l e” " š in(b,j) 
/N 
>3 V L + u ML nám ) 
n;i. n=i — 3,1 Me n=) 
l 1=0 
EV y isl jsl 
NJO 0,0 
9 iilinmer can Calcula toni Sande DisscussionSs 
a. Numerical Calculation of D 
The program for calculating L is a CRE 


n 
Alar to that used for AMSCRA with collective resolution. 


A Watfiv program of the AMSCRA with separate collision reso- 


lution is given in Appendix C. 
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D Simulation 
The Simulation program was Written mos 
and is provided in Appendix D. A simulation flow diagram is 


shown in’ igure e 








Open n slots, call random 


number generator n times 


b 







Determine i successes , 
k ccli ded i packets in b slots 





Save n=] packets 


k k = i 


J 
se 


Figure 3,11 Simulation Flow Diagram for AMSCRA 


with Separate Resolution. 
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c. Examples and Discussions 
The comparison between simulation and analysis 


of AMSCRA with separated resolution is given in Table II. 


Agreement between analysis and simulation is good. This 
verifies the analysis of the AMSCRA with separated resolu- 


tion. Some graphical results are given in Figure 3.12. 
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Figure, 322 AMSCRA with Separate Resolution 
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ANALYSIS AND SIMULATION RESULTS FOR TWO VERSIONS OF 


AMSCRA 
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Da COMPARISON BETWEEEN SEPARATE AND'COLLECTIVETRESOLUTTON 
If n < b, there is no difference between two methods. 
If n 2 b AMSCRA with separated resolution gives a better 
result than with collective resolution. Comparison between 
the two methods is given in Figures 3.13. 3.14, 3.15 and 
3.16 for various values of b. From these figures it can be 
seen that AMSCRA with separated resolution gives better 
result. These results are also provided in Table III. From 
Table III it can be seen that as n increases, the advantage 
of separate resolution over collective resolution also 


increases. 
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IV. CONCLUSION AND RECOMMENDATIONS 


The principal conclusions of this study on random cc = 
to the slotted channel are: 

l. The multislot collision resolution algorithm (MSCRA) 
is the fundamental case of this study. If we open nm slots 
after a collision involving n users, some slots are wasted 
and therefore the collision resolution interval becomes 
slightly longer: 

2. In general, the adaptive version of MSCRA, called 
AMSCRA and discussed in two different cases, this adaptive 
case gives better results than the MSCRA because we used 
smaller observation intervals in the AMSCRA. (In MSCRA, the 
observation interval always equals the number of collided 
messages, n) 

3. The first case of the AMSCRA is with collective reso- 
EII at the end of each collision resolution subperiod, 
we accumulated all collided users collectively and applied 
the algorithm to this “group. Results were better than for 
the MSCRA, as collisions were resolved in a shorter period. 

4. The second case is AMSCRA with separate resolution. 
In this case the n users involved in collision and divided 
into two subgroups. The first subgroup is resolved while 
holding the the second subgroup. Finally the second subgroup 
is resolved. Using this algorithm the best performance was 
obtained. As the performance of AMSCRA depends on the value 
of b, selecting the b=1 provided the optimum performance. 

5. The agreement between simulation and analysis is good 
for our Ewo protocols. These results could therefore be 


used to calculate throughput. 


12 


$JOB 


C*** 


a 


© 


APPENDIX A 


ANALYSIS PROGRAM FOR AMSCRA WITH COLLECTIVE RESOLUTION 


THISMEROCRAM SOLVES THE AMSCRA WITH COLLECTIVE 


RESOLUTION PROBLEM.STORAGE SHOULD BE AT LEAST 1M. 
INTEGER XX,SS,GG,QO,KK,JJ,RR,B,1,J,BB,II,Z,G,FF,C 
REAL“8 KL(O30, 30), PP(O30,030), LHAT(O30),N 
REAL*8PROB(31,31,31),SUM,MULT,SUB 
REAHS8WA(51,31,31),LHATT,L,M,P,WWV,ATZ,CR,ZS 
REAL*8 JJJ,PPP,LLL, LDENOM, LN, KKK,AAA,ZZ,TT,VV 
PRINT 10 

PRINE 19; Nime' B'. 'LN)' 

KL(1,1)=O.O 

Kez 71 )=0.0 

PRINT 99,0,0,KL(1,1) 

EPINE 99,1,0,KL(2,1) 


C ***** CALCULATION OF P(N,I) xxx 


© 


LHAT(1)=1.0 


60022... TAB, J I) COEFFICIENT *** 


© 


DO 25 BB=1,31 
DO, 26 0 zde" 

T 

NEAR I I „6541 ) 188 

ວະ Ele S pu q OR II- CT. BB) THEN 
A(BB AJA. V) <050 

ELSE IF (II.EO.JJ.AND.II.EO.BB) THEN 
FACTISFACTN(II-1) 
ABB LI) -BACH 

ELSE IF (II.EO.JJ) THEN 


TS 


FACTB-FACTN(BB-1) 

FACTBI=FACIN(BB-11) 

A(BB,JJ, II )=EAGHE Ak CTE 
ELSE 

N= 1 

ZZ=BB 

We 

Ir 

AAA=( (ZZ-VV)**(TT-VV)) 

WHILE (N.LE.(JJ-II)) DO 
AAA=AAA-A((BB-11+1),(JJ-11+1),N+1) 
N=N+1 

END WHILE 

FACTB=FACIN(BB-1) 

Fe KEBE ETI NEM, 

FACTI- FACINI E, 

FACTBI=FACTN(BB-II) 

FACTJI=FACTN(JJ-II) 

A(BB;JJ, IT)=(1./FACDE ) HA CTE 4... 


S * (1. /HACESL) ERC AZ 
END IF 
IF (((JJ-1).NE.O).AND.(BB.GE.JJ)) THEN 
CR=JJ 


ATZ- (CRI) er 
PROB (BB, JI, II)= (102222002) 
END IF 
TI າດ 
END WHILE 
26 CONTINUE 
25 CONTINUE 
C 
C*****CALCULATION OF L(N) *** 
© 
DO 61 K=2,30 
DO 60 B=1,30 
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IF (K.LE.B) THEN 
© 
C***** MESSAGE NUMBER LESS THAN OR EOUAL SLOT NUMBER*** 
C***** USE FIRST FORMULA = 
C 
SS=K 
SUM=0.0 
RR=SS-1 
ENEE EE 
MULT=PROB(SS+1,SS+1,I+1)*LHAT(SS-1) 
SUM=SUM+MULT 
5 CONTINUE 
ZS=SS 
KL(SS,B)=(1./(1.-PROB(SS+1,SS+1,1)))*(ZS+SUM) 
LHAT(SS)=KL(SS,B) 
PRINT 99,SS,B,KL(SS,B) 
ELSE 
a 
C******* MESSAGE NUMBER GREATER THAN SLOT NUMBER *** 
Ckkkkkkak USE SECOND FORMULA U AEREA AEKA 
S 
M=K 
FACTK=FACTN(K) 
P=B 
E-0.0 
poe Jar K 
FACTJ=FACTN(J) 
FACTKJ=FACTN(K-J) 
VV=FACTK/ (FACTJ*FACTKJ) 
Z=MINO(B,J) 
ES 20 1 IZ 
IF ((K-I) EO l) THEN 
ກງາກ” 
ELSE 
IHATT-ME(K-I,E) 


FS 


20 
Sch 


4O 


6O 
Si 
T5 
LO 
se 


END IF 
L=L+((A(B+1,J+1,1+1)/⁄(P**J7))*VV*((P/M)**J) 
*((1.0-(P/M))**(K-J))*LHATT) 
CONTINUE 
CONTINUE 
SUB=B+L 
LLL=O.O 
G=K+1 
DO EI 
FACTK=FACTN(K) 
FACT IN SEEN (Sa) 
EAGREJEC SE GEN(K- (J 155 
VVVEFACTK/(FACTJI"FACTJK) 
PPP=P**(J-1) 
JJJ=(PZM) FF ( Je ]) 
KO ES ey 
DLL=LLL+(A(B+1,J,1)/PPP)*VVV*JJJ*KKK 
PRINT,LLL 
CONTINUE 
LDENOM=1.0/(1.0-LLL) 
LN=LDENOM* SUB 
KE (KMEN 
PRINT 99,K,B,LN 
END IF 
CONTINUE 
CONTINUE 
FORMAT(“ *,02X,A1,12X,A1, 132mm 
FORMAT ('1') 
FORMAT(' ',13,10X,13,10X,F8.4) 
STOP 
END 


* ktp A A +“ EN NE NEN FE Ne 


*k*k FUNCTION SUBPROGRAM TO CALCULATE FACTORIAL N* 


KRERRRRKRREKRKRERRRRRKRRE RRR RAR RRR SS S SS A AS ke ka Dre 


FUNCTION FACTN(N) 
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SENTRX 


ĦAGENEL. 

IF (N.GT.1) THEN 
Bo 21-24) 

FACTN=FACTN* I 

CONTINUE 

END IF 

RETURN 

END 


7 


CIAO A O AN AO A A A O O. OSOS O 


APPEND Is 


SIMULATION FOR AMSCRA ດາ. ເໃດກ E ke ໃກາ ON 


** THIS PROGRAM WAS WRITTEN. TO MAKE EEE S 1MUB HON ໃນ? 


AMSCRA WITH COLLECTIVE RESOLUTION. 


*** VARNEJE BEN DELE IRON = x 


K 


NR 
SUM 
SUCE 
COL 


NKK 
IR 


S 
ED 
ZK 
ab 
LN 


= INPUT PARAMETER OF DISCRETE UNIFORM DISTRIBUTION 
THE lì hon E2... K GO Us SI Imi so 
PROBABILITV. K MUST BE POSITIVE IN GGUD. 

=OBSERVATION INTERVAL. 

—INPUT NUMBER OF RANDOM NUMBERS TO BEGENERATED. 

=DUMMY VARIABLE 

=SUCCESSFUL MESSAGE NUMBER 

=CURRENT COLLIDED PACKET NUMBER 

=INTERMEDIATE COLLISION RESOLUTION INTERVAL 

=DUMMY VARIABLE EQUALS B 

=DUMMY VARIABLE SHOWS SUCCESS 

=OUTPUT VECTOR OF LENGTH NR CONTAINING THE 

UNIFORMLY DISTRIBUTED INTEGERS. 

=COLLISION RESOLUTION INTERVAL IN EACH STEP 

=ENERGY DETECTOR 

=DUMMY VARIABLE 

=DUMMY VARIABLE 

=COLLISION RESOLUTION INTERVAL 


*** VARIABLE DECLARATION 


INTEGER K,NR,B, SUM, SUCC,COLL,L,Q,XX,NKK 
INTEGER IR(1),S(20000),ED(30),LAST,ZK,T 
REAL*8 LN 

DOUBLE PRECISION DSEED 

DSEED=3.DO 

Z=15000 


ze 


De Mako an 
BO a 0 
LN=O 
tO 
DĐ A = lh 
BE 
SSB 
K=Q 
GBEL-K 
SUM=K 
mel 
105 K=SUM 
IF (SUM.LE.XX) THEN 
B-SUM 
ENDIF 
NR=1 
ວສລ], K 
CALL GGUD(DSEED,K,NR, IR) 


S(I)=IR(1) 
235 CONTINUE 
DOS np 
ED(AA)=0 
301 CONTINUE 


MANS 
O ENERGY DETECTOR 
OTON OT M O So Or aa da5, 167 T7, 


* IEWIGW20721,22,28724,25.26,27,28,729, 
* 30),S(T) 
1 ED(1)=ED(1)+1 
GOTO 99 
2 ED(2)=ED(2)+1 
GOTO 99 
3 ED(3)=ED(3)+1 
GOTO 99 
4 ED(4)=ED(4)+1 


79 


10 


ija 


je? 


T3 


14 


iS 


16 


i 


18 


ປ 


20 


Z 


GOTO 99 
ED(5)=ED(5)+1 
Goro SIS 
ED(6)=ED(6)+1 
GOTO 99 
ED(7)=ED(7)+1 
GOTO 99 
ED(8)=ED(8)+1 
GOTO 99 
ED(9)=ED(9)+1 
GOTO 99 
ED(10)=ED(10)+1 
GOTO 99 
ED(11)=ED(11)+1 
"cac 

BY 0 1 
GOTO 99: 

PDS) =D) 
Comune! 99 
ED(14)=ED(14)+1 
GOTO 99 
ED(15)=ED(15)+1 
GOTO 99 
ED(16)=ED(16)+1 
GOTO 99 
ED(17)=ED(17)+1 
GOTO 99 
ED(18)=ED(18)+1 
GOTO 99 
ED(19)=ED(19)+1 
GOTO 99 
ED(20)=ED(20)+1 
GOTO 99 
EDI ໃໃ ວາ) າ 
6010 99 
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2 


Z 


24 


25 


26 


27 


28 


zo 


20 


29 


OZ 


103 


ED(22)=ED(22)+1 
GO TO 99 
ED(23)=ED(23)+1 
GOTO 99 
ED(24)=ED(24)+1 
GO TO 99 
ED(25)=ED(25)+1 
6) “ley ele 
ED(26)=ED(26)+1 
GOTO 99 
EDU ED(27)+1 
GO TO 99 
ED(28)=ED(28)+1 
GOTO 99 
ED(29)=ED(29)+1 
so Mo 99 
ED(30)=ED(30)+1 
GOTO 99 
CONTINUE 
5/6). 192 a IK 
CONTINUE 
SUM-0 
NKK=O 
D0102 T-1,B 
IE ( ED(T).EO.1) THEN 
SUCC=SUCC+1 
NKK=NKK+1 
END IF 
SUM=O - SUCC 
CONTINUE 
L=L+B 
ZK=ZK+1 
IF (SUM.EQ.O) THEN 
I0 s0) GO TO 234 
SUM-COLL-SUCC 


SU 


234 


Jill 
no 
217 


IF (SUM.LE.XX) GOTO 105 


B=XX 
GOTO 105 
ENDIF 
GOTO 105 
CONTINUE 
ENEM z 


WRITE (6,217) 0,XX,LN 

WRITE (8,217) 0,XX,LN 

CONTINUE 

CONTINUE 

FORMADI 713 5X 13. 5X Fih mD 
STOP 

END 
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$JOB 
C *** THIS PROGRAM WAS WRITTEN TO CALCULATE THE L(N) 


C 
C 
C 


C 


16 
15 
TO 


APPENDIX C 


ANALYSIS PROGRAM FOR AMSCRA WITH SEPARATE RESOLUTION 


VALUES FOR THE AMSCRA WITH SEPARATE RESOLUTION. 
STORAGE SHOULD BE AT LEAST 1M. 


MWEGER RR,B,1/J, BB IIAZACĦEEBCAUENZI,KI 
TECER KX, SS, GG, 00, KK JJ 

REAL*8 KL(030,30),PP(030,030),LHAT(030),N 
REAL*8 PROB(31,31,31),SUM,SUB,MULT,LHATT2 
REAL*8 A(31,31,31),LHATT,L,M,P,VVV,ATZ,CR,ZS 
REALF8 JJJ,PPP,LDENOM,LN,KKK,AAA,ZZ,TT,VV,LIL 
REAL*8 L2L,L12L, VBNO, VBOO 


PRINT 10 

PRINS, N', BU, LUN) 

PRINT 16,'------------------------------ ' 
FORMAT(' ',02X,A30) 

FORMAT(' ',02X,A1,12X,A1,13X,A4) 


FORMAT ('1') 
KEB: 
EEN 
PRINEN99,0,0);, KE (1, 1) 
PRINT 99,1,0,KL(2,1) 


E 222 2 CALCULATION OF HI) SA ARA AA A A K KA NAAA AK 


C 
C4 
ES 


9; 


CONTINUE 
CONTINUE 
SEET 


C*** A(B,J,I) COEFFICIENT ******* 
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DO 25 BB=1,31 
DO 26 J0 375 
La 
WEDEBEMGĦITGE.1) DO 
E (11.-CT.JJ.OR.I1.CT. BB) THEN 
Ree I e 
ELSE IF (II.EO.JJ.AND.II.EO.BB) THEN 
BAC BA LII) 
A(BB, JJ, IT) BAS 
ELSE IF (II.EQ.JJ) THEN 
FACTB=FACTN(BB-1) 
FACTBI=FACTN(BB-II) 
A(EB, JJ, TI )|=FACTE /ERETÈT 
ELSE 
N= 1 
ZZ=BB 
M 
540 
AAA=((ZZ-VV)**(TT-VV) ) 
WHILE (N.LE.(JJ-II)) DO 
AAAZAAASA (BB SI 1471) 1) 
N=N+1 
END WHILE 
FACTB=FACIN(BB-1) 
FACTI- FACH SS 
FACTI=FACIN(II-1) 
FACTBI=FACTN(BB-II) 
FACTJI=EAGHN າ 
A(BB,dJJ,11I)=(1 PC I) FACE Fin) 
*(1./FACTBI)*FACTI*AAA 
END IF 
IF (((JJ-1).NE.O).AND.(BB.GE.JJ)) THEN 
omeni 
ATZ- (CR 


84 


© 


PROB(BB,JJ,II)=(1./ATZ)*A(BB,JJ,II) 
END IF 
Mm 
END WHILE 
CONTINUE 


25 CONTINUE 


C ****CALCULATION OF L(N) *** 


© 


E 


DO 61 K=2,30 
DO 60 B=1,30 


(K.LE.B) THEN 


C***** MESSAGE NUMBER LESS THAN OR EQUAL SLOT NUMBER 
`. USESSFIRSE FORMULA 


E 


E 


SS=K 

SUM=O.O 

RR-SS-1 

POMO JRR 
MULT=PROB(SS+1,SS+1,I+1)*LHAT(SS-I) 
SUM=SUM+MULT 

CONTINUE 

ZS=SS 

KL(SS,B)=(1./(1.-PROB(SS+1,SS+1,1)))*(ZS+SUM) 

LHAT (SS)=KL(SS,B) 

PRINT 99,SS,B,KL(SS,B) 


ELSE 


E _.pO. AC GB NUMBER GREATER THAN SLOT NUMBER *** 
CX****** USE SECOND FORMULA Oki eske ke kak ke ok ake ke ak ke kk kk kok kok 


© 


M=K 
FACTK=FACTN(K) 
P=B 


85 


° 


39 


LIL-OĦO 
UL-MINO(B,K) 
BONE ME MU 


jem 
IF ((K-I).EO.1) THEN 
LHATT=1.O 
ELSE 
LHATI=KL(K-I/B) 
ENDIF 
L1L=EIL+(A(B+l,J+1,I1+1)/(P°*K)) Dp P K) 
CONTINUE 
L2L=0.0 
KÈ a 


DO 32 J=1,K1 


EAG J EAGENŲ J) 
FACTKJ=FACTN(K-J) 
VV=FACTK/ (FACTJ*FACTKJ) 
Z=MINO(B,J) 
Z1=Z+1 
DO 330n ng z1 
IF ((J-IH1).EQ.O) THEN 
EHATT=0 90 
ELSE IF (J-I+1.EO.1) THEN 
HAT O 
ELSE 
LHATTI=KÉ (JSME) 
END IF 
IF ((K-J) EO l) THEN 
LHAT O 
ELSE 
LHATT2=KL(K-J,B) 
END IF 
L2L=L2L+(A(B+1,J+1,1)/(P**J))*VV* ((P/M)**J) 
*((1.0-(P/M))**(K-J) )*(LHATT+LHATT2) 


CONTINUE 


86 


92 


6O 
Sul 
> 


SENTRY 


CONTINUE 

Ll pop hua por 
VBNO=(A(B+1,K+1,1)/(P**M))*((P/M)**K) 
VBOO=A(B+1,1,1)*((1.-(P/M))**K) 
LDENOM=1. ⁄(1. -VBNO -VBOO) 
LN=LDENOM*L12L 

KL(K,B)=LN 

PRINT 99,K,B,LN 

END IF 

CONTINUE 

CONTINUE 

FORMAT(' ',I3,10X,I3,10X,F8.4) 
STOP 

END 


* * * < * & & < & < < < k k A K A KA A * < < k < < & < < < k < k k k < < < k < < < < < < < < < 


*** FUNCTION SUBPROGRAM TO CALCULATE FACTORIAL N *** 
+ % % % ໄ % % % kk de < < # + < < + % < % % < % < % % + < < % % < < + % % < + % < % + % ແ % % % % ຫ ຫແ 
EUNCTION EACTN(N) 
BAE 26 
IF (N.GT.1) THEN 
E5 2 I=2;N 
FACTN=FACTN* I 
CONTINUE 
END IF 
RETURN 
END 


sơ 


SI 


209 EMS 

AMSC 
EE 
K = 


NR 
SUM 
SUCC 
COL 


NKK 
IR 


S 

ED 

ZK 

L 

LN 
ALAST 
CSAKAT = 
C 


A O QO GÌ (Gì OO OL OO ONO GO) O: O (200102 Q O 


DE BED ao 
MULATION FOR AMSCRA TWIETH SEPARATE TRESOEUMEON 


PROGRAM WAS WRITTEN TO MAKE THE SIMULCATITONTOE THE 
RA WITH SEPARATE RESOLUTION. 
ABLETBEEINEFFONTE A 
INPUT PAREMETER OF DISCRETE UNIFORM DISSEGNBULEION 
THE INTEGERS 1,2,...,K OCCUR Wit ESE E 
PROBABILITY. K MUST BE POSITIVE INV GUBE 
=OBSERVATION INTERVAL. 
=INPUT NUMBER OF RANDOM NUMBERS TO BEGENERATED. 
=DUMMY VARIABLE 
=SUCCESSFUL MESSAGE NUMBER 
=CURRENT COLLIDED PACKET NUMBER 
=INTERMEDIATE COLLISION RESOLUTION INTERVAL 
=DUMMY VARIABLE EQUALS B 
=DUMMY VARIABLE SHOWS SUCCESS 
=OUTPUT VECTOR OF LENGTH NR CONTAINING THE 
UNIFORMLY DISTRIBUIED Vi weeks. 
=COLLISION RESOLUTION INTERVAL IN EACH STEP 
=ENERGY DETECTOR 
=DUMMY VARIABLE 
=DUMMY VARIABLE 
=COLLISION RESOLUTION INTERVAL 
=SECOND PART OF THES COLLIDED TELSEN J 
SAVED COLLIDED .ແ ລແ, 


C *** VARIABLE DECLARATION 


E 


INTEGER K,NR,B, SUM, SUCC, COLL, L,Q, XX, NKK, FL 
INTEGER IR(1),S(20000),ED(30),LAST,ZK,T 
INTEGER ALAST(100),SAKAT(100) 

REAL*8 LN 
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DOUBLE PRECISION DSEED 
DSEED=3.DO 
Z=15000 
DO 110 XX=1,1 
DO 111 Q=2,30 
o? 
Q=4 
LN=O 
Ė-O 
DO 234 A=1,Z 
SECH 
B=XX 
SUCC=0 
K=Q 
COLL=K 
SUM=K 
ZK=1 
105 K=SUM 
IF (SUM.LE.XX) THEN 
B=SUM 
ENDIE 
NR=1 
DO 235 I=1,K 
SEENEN CDS EE E ME. RI 


EFIR 
235 CONTINUE 
DC WØO1 AA-=i, 30 
ED(AA)20 
301 CONTINUE 


POI T- 
CORO Give >, 6, /78,9,10,11,12,13 44, 15,16, 


k IME] 20,21, 22,23,24,25,26427,28,29, 
* SOM (1) 
i ED(1)=ED(1)+1 
GOTO 99 


89 


10 


dl 


1 


15 


14 


15 


16 


d 


18 


lle, 


ED(2)=ED(2)+1 
GOTO 99 
ED(3)=ED(3)+1 
GOTO 99 
ED(4)=ED(4)+1 
coto es 
ED(5)=ED(5)+1 
GOTO 99 
ED(6)=ED(6)+1 
GOTO 99 
ED(7)=ED(7)+1 
GOTO 99 
ED(8)=ED(8)+1 
GOTO 99 
ED(9)=ED(9)+1 
GOTO 99 
ED(10)=ED(10)+1 
wa 99 
ED(11)=ED(11)+1 
GO TO 99 
ED(12)=ED(12)+1 
GOTO 99 
ED(13)=ED(13)+1 
BA wa SE 
ED(14)=ED(14)+1 
GOTO 99 
ED(15)=ED(15)+1 
GOTO 99 
ED(16)=ED(16)+1 
GOTO 99 
ED(17)=ED(17)+1 
come 39 
ED(18)=ED(18)+1 
GOTO 99 
ED(19)=ED(19)+1 


JO 


20 
Zi 
22 
23 
24 
25 
R 
27 
28 
29 
30 
29 


iko) Z 


GOTO 99 
ED(20)=ED(20)+1 
GOTO 99 
ED(21)=ED(21)+1 
ejo) JE 
EB22)=ED(22)+1 
EGO! TO 99 
ED(23)=ED(23)+1 
GOTO 99 
ED(24)=ED(24)+1 
GO TO 99 
ED(25)=ED(25)+1 
se. 10 99 
ED(26)=ED(26)+1 
COTO 99 
ED(27)=ED(27)+1 
CO TO 99 
ED(28)=ED(28)+1 
GOTO 99 
ED(29)=ED(29)+1 
GO TO 99 
ED(30)=ED(30)+1 
GOTO 99 
CONTINUE 
Po 102 T=1,K 
CONTINUE 
SUM=O 
NKK=O 
POOS TIB 

IF ( ED(T).EQ.1) THEN 

SUCC=SUCC+1 


NKK=NKK+1 
ELSE 

SUM=SUM + ED(T) 
END IF 


Sa 


LOS 


106 


234 


ປ 
110 
ZAŘ 


CONTINUE 

LAST=K -NKK- SUM 

CONTINUE 

SAKAT(FL)-K-LAST 

IF (LAST.NE.O) THEN 
ALAST(EL)=LAST 
FL=FL+1 

END IF 

L=L+B 

ZKEZKt1 

IF (SUM.EO.O) THEN 
TE SUCE EO.COLL) GO m0 OE! 

SUM-ALAST(FL-1) 


EL-FL-i 
E (SUM LE. XX). COTO MOS 
B=XX 
GOTO 105 
ENDIE 
GOTO 105 
CONTINUE 
LN=L/Z 
WRITE (6,217) Q,XX,LN 
CONTINUE 
CONTINUE 
FORMAT (' ', I3,5X,I3,5X,F12.4) 
STOP 
END 
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